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In our previous work, we studied the polar distribution of the scaled energy of bottom-flavored 
hadrons from polarized top quark decays i(t) —» W + + b(—¥ XQ, using two different helicity co¬ 
ordinate systems. Basically, the energy distributions are governed by the unpolarized, polar and 
azimuthal rate functions which are related to the density matrix elements of the decay t(t) —> W + +b. 

Here we present, for the first time, the analytical expressions for the 0{a s ) radiative corrections to 
the differential azimuthal decay rates of the partonic process t(j") —> b + W + (+g) in two helicity sys¬ 
tems, which are needed to study the azimuthal distribution of the energy spectrum of the B-hadron 
produced in polarized top quark decays. Our predictions of the hadron energy distributions enable 
us to deepen our knowledge of the hadronization process and to determine the polarization states 
of top quarks. 

PACS numbers: 14.65.Ha, 13.88.+e, 14.40.Lb, 14.40.Nd 


I. INTRODUCTION 

In the Standard Model (SM), the top quark has a short 
lifetime (~ 0.5 x 10 -24 s (H) so decays rapidly and this 
short time does not allow the top quark to form the QCD 
bound states, phrased in a different language, its short 
lifetime implies that it decays before hadronization takes 
place. If it was not for the confinement of color, the 
top quark could be considered as a free particle and this 
property allows it to behave like a real particle and one 
can safely describe its decay in perturbative theory. In 
fact, at the top mass scale the strong coupling constant 
is small, a s (rrit ) ~ 0.107, so that all QCD effects involv¬ 
ing the top quark are well behaved in the perturbative 
sense. Due to the Cabibbo-Kobayashi-Maskawa (CKM) 
mixing matrix element Vtb = 0.999 [2], the decay width 
of the top quark is almost exclusively dominated by the 
two-body channel at the lowest order where a W-gauge 
boson and a bottom quark are contributed. As it is well 
known, bottom quarks produced hadronize before they 
decay (b —> JQ), therefore each Xf,-jet contains a bottom- 
flavored hadron which most of the times is a B-meson. 
The bottom hadronization is indeed one of the largest 
sources of uncertainty in the measurement of the mass 
of top quark at the CERN Large Hadron Collider (LHC) 
[3] and the Tevatron Q, as it contributes to the Monte 
Carlo systematics. The LHC is a superlative top factory, 
which allows us to carry out precision tests of the SM 
and, specifically, a precise measurement of the top quark 
properties such as its mass mt, total decay width T t and 
branching fractions. At the LHC, of particular interest 
is the distribution in the energy of meson produced in 
the top quark rest frame, so that this energy distribution 
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provides direct access to the bottom fragmentation func- 
tions (FFs). In @], we studied both the B-meson energy 
distribution produced from unpolarized top decay and 
we studied the angular distribution of the W-boson de¬ 
cay products in the decay chain t —> bW + —» Bl+vi + X. 
Since the top quark decays rapidly so that its life time 
scale is much shorter than the typical time needed for 
the QCD interactions to randomize its spin, therefore its 
full polarization content is preserved when it decays and 
passes on to its decay products. Hence, the polariza¬ 
tion of the top quark will reveal itself in the angular de¬ 
cay distribution and can be studied through the angular 
correlations between the direction of the top quark spin 
and the momenta of the decay products, IT + -boson and 
b— quark. In [f|, we studied the 0(a s ) angular distribu¬ 
tion of the scaled energy of the B-hadrons, by calculating 
the polar angular correlation in the rest frame decay of 
a polarized top quark into a stable IF + -boson and B/D- 
hadrons. We analysed this angular correlation in a helic¬ 
ity coordinate system (system 1) where the event plane, 
including the top and its decay products, is defined in 
the ( x , z) plane with the z-axes along the bottom quark 
momentum. In this system the top polarization vector 
was evaluated with respect to the direction of the bot¬ 
tom quark momentum. Basically, to define the planes 
we need to measure the momentum directions of the pt, 
and pw and the polarization direction of the top quark, 
where the evaluation of the momentum direction of pi, 
requires the use of a jet finding algorithm, whereas the 
top spin direction must be obtained from the theoretical 
input. For example, in e + e _ interactions the polarization 
degree of the top can be tuned with the help of polarized 
beams [ 7 ], so that a polarized linear e + e _ collider may be 
considered as a copious source of close to zero and close 
to 100% polarized tops. 

In Q, we analysed the 0(a s ) polar distribution of the 
B-hadron energy in a different helicity coordinate sys- 
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tem (system 2) where, as in (6| the event plane is the 
(x, z) plane but with the z-axes along the VF + -boson , so 
the polarization direction of the top quark is evaluated 
with respect to the W + momentum vector. This election 
makes the calculation so complicated. 

The azimuthal correlations between the event plane and 
the intersecting ones to this plane evaluated in two helic- 
ity systems belong to a class of polarization observables 
involving the top quark in which the leading-order (LO) 
contribution gives a zero result in the SM, so the non-zero 
contributions can either arise from higher order SM ra¬ 
diative corrections or non-SM effects [9|. Since highly po¬ 
larized top quarks will become available at hadron collid¬ 
ers through single top production processes, which occur 
at the 33% level of the tt pair production rate , it will 
then be possible to experimentally measure the azimuthal 
correlation between the (pw,Pb) and (P*,pb) planes in the 
helicity system 1 and the ( Pw,Pb ) and ( Pt,Pw ) planes in 
the helicity system 2. Here, P t stands for the polarization 
vector of the top quark and pw and pb stand for the four- 
momenta of W boson and bottom quark, respectively. To 
analyse the aforementioned azimuthal correlations in the 
polarized top rest frame, we study the azimuthal distri¬ 
bution of the scaled energy of B-hadrons at the process 
t(t) —> W + + B + X at NLO, by calculating the az¬ 
imuthal decay distribution of a polarized top quark in 
the partonic process f(t) —b b + W + /+g) in two afore¬ 
mentioned coordinate systems. For the nonperturbative 
part of the process (6, g — > B + X), from Ref. GH we ap¬ 
ply the realistic (b, g) —> B FFs obtained through a global 
fit to e+e- data from CERN LEP1 and SLAC SLC. 
Finally, we shall present and compare our numerical re¬ 
sults in both systems. These measurements will be im¬ 
portant to deepen our understanding of the nonpertur¬ 
bative aspects of B-hadrons formation and to test the 
universality and scaling violations of the B-hadron FFs. 

This paper is structured as follows. In Sec. UD we intro¬ 
duce the angular rate structure by defining the technical 
details of our calculations. In Sec. mu our analytic results 
for the 0(a s ) QCD corrections to the azimuthal distri¬ 
butions of partial decay rates are presented. In Sec. IIVI 
we shall make our predictions of energy distribution of B- 
hadrons and present our numerical analysis. In Sec. El 
our conclusions are summarized. 


II. ANGULAR STRUCTURE OF PARTIAL 
DECAY RATE 

In the current-induced t —» b transition, the dynamics 
of the process is embodied in the hadron tensor H oc 
^2 X {X\J v \t), where the SM current combina¬ 
tion is given by J / = — J//. Here, the left-chiral com¬ 

ponents of the weak current are given by oc 't/’bjpipt 
and J// oc In the transition t —>• W + + b(+g), 

the intermediate states are \X >= \b(pb) > for the Born 
term and virtual contributions and \X >= |b + g > for 


the 0{a s ) real contributions. 

The general angular distribution of the differential de¬ 
cay width dr / dx of a polarized top quark decaying into a 
jet Xb with bottom quantum numbers and a W + boson 
is expressed by the following form 

dr i (dv A dr B „ 

- = —<-b P -cos Op 

dxid cos dpd(j)p Air ( dxi dxi 

-\~P si n Q p CO s (bp 1-, (1) 

dxi I 


where the polar and azimuthal angles dp and (ftp show 
the orientation of the plane including the spin of the top 
quark relative to the event plane (see Q) and P is the 
magnitude of the top quark polarization, so P = 0 stands 
for an unpolarized top quark while P = 1 corresponds 
to 100% top quark polarization. In the notation above, 
dT a/ dx corresponds to the unpolarized differential decay 
rate, while dTp/dx and dTc/dx describe the polar and 
azimuthal correlation between the polarization of the top 
quark and its decay products, respectively. 

We shall closely follow the notation of Q, where the par- 
tonic scaled energy fraction Xi is defined as 


2 Pi ■ Pt 
m t 


( 2 ) 


As we demonstrated in Q, the finite-mt, corrections are 
rather small and thus to study the scaled energy distribu¬ 
tions of the B-meson, we employ the massless scheme or 
zero-mass variable-flavor-number (ZM-VFN) scheme (l3 | 
in the top quark rest frame, where the zero mass parton 
approximation is also applied to the bottom quark. The 
non-zero value of the b-quark mass only enter through 
the initial condition of the nonperturbative FFs. Non¬ 
perturbative FFs are describing the hadronization pro¬ 
cesses (6, g) —>■ Xb and are subject to Dokshitzer-Gribov- 
Lipatov-Alteralli-Parisi (DGLAP) evolution [L3j. 

By the zero mass approximation, one has 0 < Xi < 1 — oj 
where w = rn/^/m/. Throughout this manuscript, we 
apply the normalized partonic energy fraction as 


2 Ej 

m t (l - u) ’ 


(* = b,g) 


( 3 ) 


where Ei refers to the energy of outgoing partons (bottom 
or gluon) and 0 < Xi < 1. 

In our previous works, the NLO radiative corrections 
to the unpolarized differential rate dT a/ dxi Q and the 
polar differential rates dTp/dxi j(3, [§J have been stud¬ 
ied in two possible helicity systems, extensively. In the 
present work, we study the radiative corrections to the 
azimuthal correlation function dTc/dxi in both helicity 
systems, which have not been done before. Finally, at 
the hadron level we shall compare our predictions for the 
energy distribution of B-mesons in two coordinate sys¬ 
tems 1 and 2, considering all contributions. 
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III. ANALYTIC RESULTS FOR AZIMUTHAL 
DECAY DISTRIBUTIONS 

In the rest frame of a top quark decaying into a W + bo¬ 
son, a b-quark and a gluon, the final state particles define 
an event plane. Relative to this plane one can then define 
the spin direction of the polarized top quark. There are 
two different choices of possible coordinate systems rela¬ 
tive to the event plane where one differentiates between 
lielicity systems according to the orientation of the z-axis 
(these systems compared in Q). Four-momenta of the fa- 
quark and the W + boson in these two various coordinate 
systems are defined as 

System 1 : p b \\z \{pw)x > 0 

System 2 : pw\\z ]{Pb)x>0 (4) 

Indeed, in the system 1 the three-momentum of the fa- 
quark points into the direction of the positive z-axis and 
in the system 2, the three-momentum of the W + boson 
is defined along this axis. 

In the following, we explain the technical details of our 
calculation for the NLO radiative corrections to the tree- 
level decay rate of t(t) —> b + W. 


A. Born term results 

In the SM, the polarized top decay rate is dominated 
by the decay process t(p t ) —> b{p b ) + W + (pw) at the 
Born level. In the rest frame of the top quark, the four- 
momentum of the top quark is set to pt = (mp, 0) and 
the polarization four-vector of the top quark is set as 
St = P(0; sindp cos (j)p, sin 9p sin 0p, cos 9p), where P is 
the top polarization degree (0 < P < 1). Considering 
the coordinate system 1, where the three-momentum of 
the b-quark points into the positive z-axis, we set the 
four-momentum of the b-quark as p b = E b ( 1; 0,0,1) and 
in the system 2, it is p b = E b (l; 0,0, — 1) where the three- 
momentum of the W + boson is defined along the posi¬ 
tive z-axis. Note that by applying the ZM-VFN scheme 
we put the b-quark mass to zero throughout this paper. 
Thus, the Born term helicity structure of partial rates, 
reads 


d 2 r(°> 

d cos 9pdcj)p 


i-jr^TPr^ cos 9p 
+Pr^ 0) sin 6» P cos ^pl, (5) 


where the sign stands for the helicity system 1 and 
the sign ,, +" is for the second one, and 


where, 9\y is the weak mixing angle and a is the tiny 
structure constant. These results are in complete agree¬ 
ment with the expressions in 0,0. 

As it is seen, the Born term contribution to Tc is zero. 
We point out that the vanishing of this azimuthal cor¬ 
relation is a consequence of the left-chiral (V-A)(V-A) 
nature of the current-current interaction in the SM. An¬ 
other example of a LO zero polarization observable is 
given in 0. There, we showed that the decay of a top 
quark into a polarized transverse-plus W boson and a 
(massless) bottom quark leads to the contribution zero 
for the top decay rate into the transverse-plus W boson 
at the Born term level due to the left-chiral (V-A) cou¬ 
pling structure of the SM. However, if one takes a massive 
b-quark in the calculation, this contribution is no longer 
zero but the LO result obtained for Tc does not depend 
on the mass of the bottom quark. 


B. QCD NLO contribution to the azimuthal 
differential decay rate dT c /dxi 

Generally, the required ingredients for the NLO per¬ 
turbative calculation are the virtual one-loop contribu¬ 
tions and the tree-graph (real emission) contributions. 
Since, at LO the relevant scalar products are pt ■ St = 0, 
Pt ■ p b = mtE b and p b ■ St = pPE b cosOp (in both he¬ 
licity systems) then the virtual one-loop corrections are 
contributed in the unpolarized rate (T^) and the polar 
correlation function (Tb), which have been studied ex¬ 
tensively before 0, while the azimuthal one (Tc) does 
not have any contribution from the virtual corrections. 

The QCD NLO contribution results from the square 
of the real gluon emission graphs. By working in the 
massless scheme where m b = 0, for the corresponding 
real amplitude squared one has 

| M real|2 = _ 7 T 2 C F Oa B + P^Pw x 

4 sin 2 9 W V m^y 

f Pi | P 2 _2 P 3 \ 

1 (PfPg) 2 ( Pb'Pg ) 2 (PfPg)(pb-Pg) )’ 

(7) 

where Cf = 4/3 stands for the color factor, and 

Pi = - 7 5 )(m t + - tf g )^{rft + m t ) x 

(1 +75?t)7/3(TOt +#t -j4)(1 + 75)7!/], 

P 2 = Trlrfb^^b +rf g )lp0--lb){m t +$)(! +75$) x 

(1 + 75 ) 7^(24 +]/ghp\, 

P 3 = Tr[pt b 7 M (i - 7 5 )(m t +y4 - j/ g )7 P{m t +rft) x 

( 1 + 75 $) ( 1 + 75)7 + tfghp\- ( 8 ) 


,(0) 

mtOL 

A 

16w sin 2 9w 

40) 

mt®. 

B 

16w sin 2 9w 

40 ) 

= 0 

c 


(1 + 2w)(l — w) 2 , 
(1 — 2w)(l — w) 2 , 


( 6 ) 


In general, to regulate the gluon IR singularities we work 
in a D-dimensions approach, where the differential decay 
rate for the real contribution is given by 


2(4 -D) 

dr = tZE -1 M 1 

2 mt 


2 dR 3 (p t ,pb,p gi pw), 


(9) 
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where, fi is an arbitrary reference mass and the 3-body 
phase space element dR 3 reads 


d D ^bd 0 1 pw-- d D 1 P g , 0 _,3-2D^/„ 

op or or (~ /l ) \Pt 

ZiU/]} —j Jl 2 j g 


Y pf)- 

g,b,W 

( 10 ) 


of the detected meson so that this contribution decreases 
the size of decay rate at the threshold energy. dT a/ dx g is 
the same in both coordinate systems and can be found in 
[ 5 }, and the analytical expression for the dTip/dxg in the 
helicity system 1 is presented in Q and the dr 2 b/ dx g is 
given in [8j, where x g is defined in ]3]). In the coordinate 
system 1, the azimuthal differential width dT/dx g reads 


Here, d D 1 |p| = \p\ D 2 d\p\dQ where the angular integral 
in D-dimensions will have to be written as 


dfl 

d(f>pd cos dp 


~ wD- 3 ' ( sin ^) fl 4 (sin <j)p) D 4 . 

1 \ 2 / 

(ii) 


Considering the general form of the angular decay distri¬ 
bution dT]), as we showed in [§f the unpolarized differential 
decay rate dT a/ dx is independent of the applied helicity 
system but the polar distribution of decay width dr p/dx 
depends on the various choices of possible coordinate sys¬ 
tems, but all final results are free of IR singularities. In 
the following we will concentrate on the differential az¬ 
imuthal correlation function dTc/dx , considering both 
helicity coordinate systems. 

In the system 1, the relevant scalar products are 


p g ■ St = —PE g (sm9 g bsindp cosc/p + cosdgbcosdp), 
Pg ' Pb — E g Eb{ 1 COS 0 g b ), 

Pb-St = —PE b cos 6p, (12) 


where 0 g b is the polar angle between the gluon and the 
bottom quark momenta in the event plane, so cos 0 g b = 
( to ? - to ^ - 2 m t (E b + E g ) + 2E b Eg)/(2E b E g ). To cal¬ 
culate the dTc/dxb, in J3]) we fix the momentum of the 
b-quark and integrate over the energy of the gluon, which 
ranges as nritS{ 1 — x b ) < E g < mtS( 1 — x;,)/(l — 2Sx b ), 
where 5 = (1 — w)/2. Therefore, one has 


dTic 

dx b 


n (0) Oi s CF 

B 2(1 — w)(l — 2ui) 
4(1 + u) 


+ 11 - 


4 (w — l)x b — 2uj 2 — 5 u> 

1 


x b 


(1 - Xb{l - w))2 L 


- 16 


2w + 10w 2 + 4 ( l+UJ ^ + 2(1 - w ) 2 (u; - 5 )x 2 b + 
x b 


dT\c 
dx a 


T 4(1 - “ )x » - 8 + 6 “ + 


2(1 + u>) 1 — 2w 


1 


X a X ‘ 4(1 -Xg(l -w))3 L 


-7(1 - w) 3 x 3 - 2(1 - w) 2 (5w - 19)x 2 - 


(9c ij + 19 u) — 93w + 65)x g — 


2 ^-'-19)x 2 
4(1 - 2w) 


3 t 2 

x g 


2(2w 3 - 13w 2 -5w + 20) + 


2(6w 2 - 13w - 1) 


(14) 


In the helicity coordinate system 2, the relevant scalar 
products are 


Pg ■ s t 
Pb ■ s t 


Pg-Pw 

Pb-Pw 
Pw ■ St 


PE g (sin0 g w sin dp cos tj>p — 
PEb (sin Obw sin 6p cos <j)p — 


m t (E w + Eg) - 


L w 


m t {E w + E b ) - 


TO? + TO 


W 


-P\pw\ cos Op, 


cos 0 g w cos Op), 
cos Obw cos Op), 


(15) 


and Pt • s t = 0. In the system 2, Obw is the polar angle 
between the b-quark momentum and the W boson (z- 
axis) and 0 g \y is the angle between the gluon and the W 
boson, whereas cos 0 g w = (to? + — 2m t (Ew + E g ) + 

2E w Eg)/(2E g p w ) and cos 6W = (to? Tto^ - 2m t (E b + 
Ew) + ZEbEw)/ (2 E b pw) with p w = — m 

As before, to calculate the azimuthal differential rate 
dTc/dx b , in ([21) we fix the momentum of the b-quark 
but we integrate over the energy of the W boson, which 
ranges as m t (uj + [1 — 2Sx b ] 2 )/ (2(1 — 2 Sx b )) < E g < 
nit (1 — S). Therefore, in the coordinate system 2 the 
azimuthal differential width dTc/dxb is expressed as 


(1 - u) 3 xl + (7w 3 - 3w 2 - 25 to + 21)x b \ j, 

(13) 

where x b is defined in ([3]). This result can be compared 
against known results presented in 0 after integrating 
over x b (0 < x b < 1). 

Since the observed mesons in top quark decays can be 
also produced through a fragmenting real gluon, there¬ 
fore, to obtain the most accurate energy distribution of 
the B meson one has to add the contribution of gluon 
fragmentation to the b-quark one. In Q, it is shown that 
the gluon contribution can be important at a low energy 


dT 2 c 


= r 


(O) 


o. s Cf 


dxb B 2(2 — (1 — bo)xb) 2 (l — w)(l — 2ui) 

\x\{ 1 - w) 2 (2w 2 - 3w - 3) + 4w 3 + 22w 2 - 16w 
—2 + 2xf,(l — w)(—6ca 2 + 7(jJ + 3) — 

o IU + (1 - Xb(l - u)) 2 r_ 2/1 ^2/0 ^ , 

2\l - 1 — w) (2cu — 1 j + 


1 + CJ 

2xjj{i0 — -|- 4c j — 1) -|- (1 -|- cj)(2cj^ -|- 9cj — 1) 


+8w(2 - Xb(l - w)) 2 y/l - Xb( 1 - w) j, 


(16) 































5 


and for the gluon one, we have 

dT 2 c _ r (o) ol s Cf 
dx n B 


(1 - w)(l - 2 u) { 2x g ^/l - x g { 1 - w) - 
4 1 - x g) ~ (! - 2u 2 )(l - x g ) 

( /^ + (! ~ x g (l - uj)) 2 


2w — 1 — 


i+„ -v 1 -^ 1 -' 

(1 + w) 2 (l - 2uj) 2(1 + 2w 2 ) 




(17) 


IV. NUMERICAL ANALYSIS 


By determining the differential decay rates in the par- 
ton level, in the first step we turn to our numerical pre¬ 
dictions of the unpolarized and polarized decay rates by 
integrating dT/dxb over Xb(0 < Xb < 1), while the strong 
coupling constant is evolved from a s [mz ) = 0.1184 to 
a s( TO t) = 0.1070. By combining our results for the dif¬ 
ferential azimuthal correlation functions (Eqs. (1131) and 
(llpp i with the results obtained for the unpolarized rate 
Q and the polar correlation rate in the system 1 Q, one 
has 


dr" LO 

dtfiipdcosdip 


(1 -0.08542) 

-l4 0) (l - 0.1303)Pcos6»ip 
—^ (0.0980)P sin 9\p cos 01P | 

pNLO , 

= -4—^ 1 - 0.3777P cos 9ip 
47t L 

—O.O426Psin0ip cost)>ip|, (18) 


and by considering the polar correlation one in the sys¬ 
tem 2 [lj, one has 


dr* LO 

d<f> 2 pdcos 9 2 p 


|r^ 0 ) (l- 0.08542) 
+T ( °\l - 0.2814)Pcos 9 2 p 


+r^ 0) (0.01446)Psin 9 2P cos <t> 2P j 

t-NLO r 

= -4—U +0.3121 cos0ip 
47r l 

+O.OO63Psin0ip cos^ipj, (19) 

where = 1.4705 and ^ = 0.5841 if one sets mw = 
80.399 GeV, m t = 172.98 GeV and sin 2 9 W = 0.2312 03- 
As it is seen, the azimuthal correlation generated by the 
radiative corrections is quite small, especially in the sec¬ 
ond coordinate system. We can assert that, if top quark 
decays reveal a violation of the SM (V-A) current struc¬ 
ture in the azimuthal correlation function which exceeds 



Figure 1: xb distribution of dT Nlj0 /d,XB in the helicity sys¬ 
tem 1, considering the unpolarized (dotted line), the polar 
(solid line) and the azimuthal (dashed line) contributions. 



Figure 2: As in Fig [I] but in the helicity system 2. 


the 5% level at the system 1 and the 1% level at the 
system 2, the violation must have a non-SM origin. 

In the last step, we present our phenomenological re¬ 
sults for the energy spectrum dr /dxp of the B meson, 
where we define the normalized energy fraction of the 
B meson as Xp = 2Eg/(m t (l — w)) 13]). According 
to the the factorization theorem of the QCD-improved 
parton model 03, the energy distribution of B meson 
can be expressed as the convolution of the parton-level 
spectrum dr /dxi with the nonperturbative fragmenta¬ 
tion functions D b (z,[1f) as 


dr 

dxp 


= E 


i=b,g 


~r~ ( 9 -r, 9 -f)®D b (— ,/rp). (20) 

ax j x j 


The integral convolution is defined as (/ <g) g)(x) = 
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fx dxf{z)g{x/z). In (1201) . ji y and ji B arc tlic factor¬ 
ization and the renormalization scales, respectively, that 
the scale g, B is associated with the renormalization of the 
strong coupling constant and a choice often made consists 
of setting = hf- As in our previous works, we adopt 
the convention = fJ-F = m fi . 

In (l20l) . Di(z,HF) is the nonperturbative FF describing 
the transition ( b,g ) —>• B which is process independent. 
Several models are proposed to describe the nonperturba¬ 
tive transition from a parton into a hadron state. Here, 
following Ref. El we employ the B meson FF deter¬ 
mined at NLO in the ZM-VFN scheme and obtained 
by fitting the experimental data from the ALEPH and 
OPAL collaborations at CERN LEP1 and by SLD at 
SLAC SLC. Authors in El have parametrized the z dis¬ 
tribution of the b —i B FF at the initial scale po = rnb as 
D^{z,g, o) = Nz a (l — z)P (power model), while the gluon 
FF is set to zero at the starting scale and is evolved to 
higher scales using the DGLAP equations El. Their 
results for the fit parameters at the initial scale are 
N = 4684.1, a = 16.87 and /3 = 2.628 with 3? = 1.495. 
Following Ref. [l6|, as numerical input values we take 
m w = 80.339 GeV, m b = 4.78 GeV, m B = 5.279 GeV 
and the typical QCD scale A^L = 231 MeV. Note that, 
in the ZM-VFN scheme the b-quark mass only enter 
through the initial condition of the FF. 

To study the x B scaled energy distributions of B mesons, 
we consider the quantity dr(t('|') —> B + X)/dx B in the 
two helicity coordinate systems. In (Hj, [f|, we showed 
that the g —i B contribution into the NLO energy spec¬ 
trum of the B-meson is negative and appreciable only in 
the \ow-x B region and for higher values of x B the NLO 
result is practically exhausted by the b —^ B contribu¬ 
tion. The contribution of the gluon is calculated to see 
where it contributes to dr/ dx B and can not be discrim¬ 
inated in the meson spectrum as an experimental quan¬ 
tity. In the scaled energy of mesons, all contributions 
including the bottom quark, gluon and light quarks con¬ 
tribute. In Fig. [TJ the a;B-spectrum of the B meson pro¬ 
duced through the unpolarized top quark decay (dotted 
line) is shown. The polar (solid line) and the azimuthal 


(dashed line) contributions in the helicity system 1 are 
also studied. As we explained, the azimuthal correlation 
is prohibited at LO, which explains the smallness of the 
corresponding result. Note that the threshold occurs at 
x B > 2m B /(m t (l — to)) = 0.078. In Fig. [2J the same 
predictions are shown in the helicity system 2. 


V. CONCLUSIONS 

To study the 0(a s ) spin-dependent energy spectrum 
of hadrons produced from polarized top quark decays, 
one needs to know the NLO radiative corrections to the 
angular differential decay rates of the process f(f) —> 
b + W + . In our previous works, the unpolarized decay 
rate (<£T a/ dxi) and the polar correlation one ( dY B /dxi ) 
were calculated at the parton-level in two different he¬ 
licity coordinate systems. These various helicity systems 
provide independent probes of the polarized top quark 
decay dynamics. 

Here, by considering two helicity systems we have calcu¬ 
lated the 0(a s ) corrections to the differential azimuthal 
correlation function ( dY c /dXi ) which vanishes at the 
Born term level. These quantities are required to calcu¬ 
late the x B distribution dVc/dx B of t —>• B + X. Com¬ 
paring future measurements of the polarized and unpo¬ 
larized partial widths dT/dx B at the LHC with our NLO 
predictions, one will be able to test the universality and 
scaling violations of the B meson FFs. These measure¬ 
ments will finally be the primary source of information 
on the B meson FFs and the azimuthal correlation func¬ 
tion dTc/dx B can also constrain the g —>• B and b —> B 
FFs even further. 

We also found that the 0(a s ) radiative corrections to 
the azimuthal correlation observable T c are so small, es¬ 
pecially in the helicity system 2. Specifically, it is safe 
to say that, if top quark decays reveal a violation of the 
SM (V-A) current structure in the azimuthal correlation 
function which exceeds the 5% level at the system 1, the 
violation must have a non-SM origin. 
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